All-optical reconstruction of atomic ground-state population 



P. London, 1 0. Firstenberg, 1 M. Shuker, 1 and A. Ron 1 
1 Department of Physics, Technion- Israel Institute of Technology, Haifa 32000, Israel 

The population distribution within the ground-state of an atomic ensemble is of large significance 
in a variety of quantum optics processes. We present a method to reconstruct the detailed population 
distribution from a set of absorption measurements with various frequencies and polarizations, by 
utilizing the differences between the dipole matrix elements of the probed transitions. The technique 
is experimentally implemented on a thermal rubidium vapor, demonstrating a population-based 
analysis in two optical pumping examples. The results are used to verify and calibrate an elaborated 
numerical model, and the limitations of the reconstruction scheme which result from the symmetry 
properties of the dipole matrix elements are discussed. 



I. INTRODUCTION 

In the past decades, alkali-metal atoms have become a 
central ingredient in quantum optics applications, mostly 
due to the simplicity of their electronic internal struc- 
ture. In many atom-photon phenomena and particularly 
in multi-photon processes, a key element is the distri- 
bution of populations among the various states within 
the alkali ground-level manifold. As an example, con- 
sider the implementation of a frequency reference based 
on coherent population trapping (CPT), in which two 
resonant fields are used to create and interrogate the co- 
herence between two specific ground-level states 0, 0- 
Dilution of the CPT states due to optical pumping and 
the redistribution of population towards highly polarized 
states pose a major limit on the CPT contrast [||. In 
addition, schemes for narrow CPT resonances have been 
proposed in which spin-exchange decay is suppressed by 
the use of a proper population distribution For such 
realizations of CPT and similar physical arrangements, 
the detailed electronic population distribution is of great 
importance. 

Population reconstruction is a measurement of the pro- 
jection of an unknown density matrix along a certain 
quantization axis (diagonal elements). If similar mea- 
surements are carried out from sufficient directions, a to- 
mographic reconstruction of the full density matrix (off- 
diagonal elements) is possible [5|. Various reconstruction 
schemes have been demonstrated in recent years 043; 
offering diverse tools for the characterization of quan- 
tum optics processes and for the verification or calibra- 
tion of theoretical models for atom-photon interactions. 
In many cases, the measurements themselves or the ro- 
tations of the measured axes perturb or even erase the 
measured state. In other cases, the requirements of the 
probing protocol contrast with the operating conditions 
of the measured system, e.g. an application of a mag- 
netic field in a magnetically-sensitive arrangement [13j . 
A continuous monitoring of the state dynamics requires 
an interrogation method that does not substantially dam- 
age it, which is possible for non-radiating systems us- 
ing weak all-optical probing. Recently, all-optical mea- 
surements of dark-resonance spectra have been used to 
extract the ground-state populations of cold rubidium 



atoms jlOj. However, since dark resonances are two- 
photon transitions requiring non-linear pumping of a cer- 
tain dark state, they cannot be carried out with very 
weak probe beams. Here we propose an all-optical lin- 
ear interrogation method with essentially no lower bound 
for the probe intensity. Having a negligible effect on the 
ensemble, the procedure can be carried out side by side 
with other processes. 

In [l4j, re-population processes within the ground- 
state manifold of a thermal rubidium vapor were inves- 
tigated using a series of weak light pulses of single fre- 
quency and polarization. The time-dependent absorption 
was analyzed to estimate the ground-state decay and de- 
coherence rates. The scheme presented here generalizes 
this method for multiple frequencies and polarizations, 
in order to characterize the detailed sub-level population 
distribution within the ground manifold. 

The populations of the different angular-mumentum 
states within the ground level of an atomic ensemble 
can be readily measured in a spectroscopic measurement, 
only when the Zeeman and the hyperfine splittings are 
both larger than the spectral widths of the transitions. 
In a room-temperature atomic vapor, the optical Doppler 
broadening is on the order of hundreds of MHz, and the 
collisional (pressure) broadening in a typical case of a 
buffer-gas environment may even exceed that. Thus, the 
magnetic lines are resolved only for magnetic fields much 
larger than 100 Gauss, which are not suitable in many 
setups, such as storage-of-light 
erences [1, EH, magnetometry 
18]. In these situations the angular-mumetum states 
are degenerate with respect to the measurement. We 
show that this effective degeneracy may be lifted utilizing 
the differences in the transition strengths for each state 
and each polarization. Therefore in contrast to previous 
schemes, the reconstruction can be carried out even when 
the contribution of any single state in each measurement 
could not be distinguished, and thus no application of a 
significant external magnetic field is required. The ba- 
sic concepts of the reconstruction method are outlined 
in Sec. II, and the experimental setup of the probing 
system is described in Sec. III. The results for various 
optical-pumping experiments are presented in Sec. IV. 
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II. THEORY OF POPULATION 
RECONSTRUCTION 

We study optical atomic transitions of the form 
\F, m F ) —> \F',rriF>) , where F and rap are respectively 
the hyperfine and Zeeman numbers within the ground- 
level manifold, F' and rriF' denote the excited mani- 
fold. For the sake of simplicity, we consider here only 
circularly-polarized light, with mpi = m F ±l correspond- 
ing to the o~ polarizations. The absorption of a weak 
probe field due to each atomic transition depends on the 
inherent strength of the transition, determined by the 
electric dipole moment, and on the frequency detuning 
from the resonance. The absorption coefficient is thus 
written as the sum over the contributions of all possible 
transitions, 



M= £ c± mF) HxP |Fil 



(i) 



\F,m F ) 



where 



C \F,rn F )(u)=n — Y J K FmF {u)(HF,m F ) ( 2 ) 



is the excitation spectrum of a specific ground-level state 
(\F, nip)), P\F,m F ) is its population, and n is the atomic 
density. The dipole matrix element is expressed in terms 
of the Clebsch-Gordan coefficients as |19| 



F',m F ±l 
^F,m F 
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(3) 



where /^o factors out the coefficients common to all tran- 
sitions, J and J' are the total electronic angular mo- 
menta of the ground and excited levels, respectively, and 
I is the nuclear angular momentum. The frequency de- 
pendence of the spectrum is introduced via the Voigt 
functions, Kp^ F±1 (uj), centered around the resonance 
frequencies of \F,mp) — > \F' ,mpr) and accounting for 
both the homogeneous and the inhomogeneous broaden- 
ing components. Sampling the spectra of the er + and o~ 
polarizations at TV different frequencies transforms Eq. 
(fTJ) into a matrix form, 



a (2iVxl) — C(2AfxAf) ' P(Mxl)) 



(4) 



where each row corresponds to a single measurement, a 
is the vector of measured absorption coefficients, P is the 
ground-state populations, and C is the coupling matrix 
obtained from Eq. ([2]). The population reconstruction 
is thus possible via a matrix inversion, providing that 
at least M linearily-indpendent measurements are made. 
A least-squares fit to the data with 2N ^> M measure- 
ments would increase the accuracy and robustness of the 
reconstruction, and is more adequate in the presence of 
real-life measurement noise and systematic errors. Such 
a fit to Eq.Q is given by 

P = C+a, (5) 



where C + = (C T C)~ 1 C T is the Moore-Penrose pseu- 
doinverse [2(J. Finally, if the atomic density, n , is not 
known a priori within the desired accuracy, it should be 
determined by the constrain £-P|F,mp) = 1< In practice, 
this can be achieved directly from the measured spectra 
using the relation 

Jduj[a + (uj) + a~(w)] = Qtiq^^-, (6) 

where the identity ~£f> I(m + ) 2 + (/O 2 ] = Q (Q = 2/3 for 
rubidium-87 Dl transitions) and J duK{ui) = 1 are used. 

In order to demonstrate the strengths and the limita- 
tions of our scheme, we hereafter concentrate on the Dl 
transition of rubidium-87. The ground and upper man- 
ifolds of this transition consist of two hyperfine levels, 
F = 1,2, which in turn are Zeeman splitted, as depicted 
in FigJTJj. When the Zeeman and the hyperfine splittings 
are larger than the spectral line- width, it is easy to show 
that C T C is nearly diagonal and trivially invertible, and 
the populations of the eight ground- level states are imme- 
diately obtained from Eq.([5]). In the opposite limit where 
both splittings are small, and the spectral lines are indis- 
tinguishable, two of the eigenvalues of C T C in Eq.© are 
close to zero due to the Clebsch-Gordan properties [27l |. 
For this limiting case, our reconstruction scheme becomes 
sensitive to noise and errors and thus not applicable. 

Nevertheless in the most common intermediate regime, 
the Zeeman sublevels are degenerate but the ground and 
excited hyperfine sublevles are spectrally resolved. For a 
thermal rubidium vapor, the ground-level hyperfine split- 
ting is much larger than the spectral line- width (6.8 GHz 
compared with ^500 MHz), and the absorption spec- 
tra of the F = 1 and the F = 2 manifolds are com- 
pletely separable. The spectrum of each F manifold, for 
a given probe polarization q, provides two meaningful 

data points, £p ~ 1,q and £p ~ 2 ' q , which are the sum of all 
(Zeeman) transitions between F and F'. A fitting pro- 
cess a q (uj) — J2f f'=i 2^f ' q ^F where Kp (u) are 
the degenerate Voigt profiles, gives £ F ' y with high ac- 

F' a 

curacy. With £ F ' H known, one is left with the following 
reduced equation, 



tF' tl ( ,,F' ,m F ±q,2 r> 

SF — l^\F,m F )\t l F,m F ) ^\F,m f 



(7) 



for each ground-state manifold. For example, using cir- 
cular polarizations , a rank-3 matrix, 
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is obtained for the F — 1 manifold, with a solution 
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FIG. 1: (color online). The setup of experiment I: a. 
VCSEL — vertical-cavity surface-emitting laser diode; PBS — 
polarizing beam splitter; PLZ — GLAN polarizer; MS — 
magnetic shield, EC — electric coil; ECDL — external-cavity 
diode laser, b. Energy levels in experiment I, and exciting 
pump fields, showing resonant transitions in red solid lines 
and non-resonant transition in blue dashed lines. The end- 
state \F = 2, tuf = +2} is circled, c. Energy levels in experi- 
ment II and exciting pump fields. 



III. EXPERIMENTAL SETUP 

The experimental arrangement, with both the pump- 
ing and the probing setups, is depicted in Fig. [Tk, We 
drive the atomic ensemble to a specific atomic state us- 
ing a pump beam, as detailed in sections IIV Al and IIVBI 
Simultaneously with the pumping process, we employ 
a continuous weak light field for probing the medium 
(Fig. QJ,). For the probe beam, we use a vertical-cavity 
surface-emitting laser diode (VCSEL), and its frequency 
is scanned by ramping the current. The probe beam, typ- 
ically of ~lmm diameter and ^50nW power, traverses a 
Glan-Laser polarizer followed by a quarter- wave plate, to 
produce a o~ + polarization. A half-wave plate, mounted 
on a flip-flop, is used to convert the probe polarization to 
(j~ in alternating measurements. The probe then passes 
through the main vapor cell and is measured by a photo- 
diode. The frequency scans are of ~12 GHz and carried 
out at a repetition rate which is sufficiently small in or- 
der to probe the stady-state distribution (12.5 Hz). Each 
scan consists of 2500 data points (See N in Eq. @). 
The population reconstruction scheme studied here re- 
quires high accuracy determination of the frequency of 
the probe. To avoid errors caused by frequency drifts of 
the VCSEL, and nonlineartities in the frequency scan, a 
frequency calibration system was utilized. An accurate 
frequency scale was acquired by a saturated absorption 
spectroscopy setup [2l[ in a reference vapor cell. Non- 
linearities in the frequency scan were corrected using fre- 
quency markers produced by a Fabry- Perot cavity (22j. 
After the frequency calibration process, the absorption 
spectra were introduced into Eq. ([5]) and regressed with 
a numerical code. 



For the F = 2 manifold, the matrix obtained is also of 
rank 3, whereas the manifold contains 5 sublevels, thus 
preventing the inversion. The addition of other excita- 
tions, either to different levels (e.g., the D2 transitions) or 
with a 7r-plorization, are ineffective in resolving this case. 
This follows from the symmetry properties of the dipole 
matrix elements, namely the absorption of an isotropic 
light field [l^| and the spherical behavior of the dipole 
moment, 



F'.m F - 
^F,m F 



F',(-m F )+q 
^F,(-rn F ) ' 



(10) 



enabling no more than three independent relations. In 
conclusion, the detailed populations of the F = 1 mani- 
fold, i.e. Pm-u, -P|i,q)> and Pn+i), can be recovered, but 
the F = 2 populations can not. The total population in 
the F = 2 manifold is determined, and in specific cases, 
as we show in what follows, the total F = 2 population 
can be related to that of the maximally-polarized states 
|2, ±2), in effect reconstructing the complete ground-level 
population distribution. 



IV. EXPERIMENTAL RESULTS 

A. Experiment I: Optical pumping with a + 
polarized light 

In the first experiment, we pump the atomic ensemble 
with an external-cavity diode laser (ECDL), stabilized to 
the F = 2 ->■ F' = 2 transition of the Dl manifold (Fig. 
[]})). As illustrated in Fig. QJ,, the pumping laser passes 
through a polarizing beam splitter (PBS) followed by a 
quarter-wave plate to produce the circular polarization. 
The beam is then directed into the 5-cm long main va- 
por cell, containing isotopically pure 87 Rb and 10 Torrs 
of neon buffer gas. The cell is at room temperature, pro- 
viding a vapor density of ~10 10 /cc. The cell is located 
within a three-layered magnetic shield, and an electric 
coil produces a small axial magnetic field (B z — 450mG) 
to set the quantization axis. 

We have performed population-distribution measure- 
ments with this arrangement, as described in Sec. IIII1 for 
various pump powers. Figure [5] depicts probe absorption 
curves for three typical pump powers. When no pump 
is applied (panel b in Fig. [5]) the medium is expected 
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157% FIG. 3: The reconstructed populations versus five decads of 
the pump intensity in experiment I (symbols). The numerical 
model (lines) agrees well with the measurements. 



FIG. 2: (color online) Measured signals in experiment I: a. 
FP signal (blue thin lines) and SAS signal (red thick line) are 
used for callibrating the frequency axis. (b,c,d) Absorption 
spectra measured with a + polarization (solid lines) and a~ 
polarization (dashed lines) with various pump powers, show- 
ing the thermal equilbrium (b), depopulated F = 2 (c), and 
vacancy of the mp < 1 states in F = 2 (d). On the right are 
the population distribution of the F = 1 manifold and the 
total population of F — 2 manifold. 



to be in thermal equilibrium, in which the population is 
equally distributed among the eight levels of the ground- 
state. Indeed, the measured F — 1 populations are equal 
up to the measurment error. As the pump power is in- 
creased, the absorption from the F — 2 manifold is re- 
duced due to the depopulation of the levels within the 
F = 2 manifold. The a + and er~ absorption curves dif- 
fer as a consequence of the increasing polarization of the 
medium. For pump intensities larger than 10/iW/mm 2 , 
the a + absorption from the F = 2 manifold nullifies, 
showing that the \F = 2,rriF = —2, —1,0, +1) levels are 
vacant. 

The populations of the three states within the ground- 
state F = 1 manifold are reconstructed using the pro- 
cedure described in Sec. HU The best-fitted values of 
the homogeneous and inhomogeneous broadenings of the 
Voigt lineshape, T =103 MHz and a =202 MHz (FWHM 
^600 MHz), are as expected for our room-temperature 
buffer-gas system. Although the population distribution 
of the F = 2 manifold is not formally solved, a least- 
squares minimization process assuming non-negative val- 
ues can be used to estimate the detailed populations 
within the F — 2 manifold. At high pump powers, 
in which the population of F — 2 is associated with 
\F = 2, mp = +2) only, this estimation becomes reliable 



and the populations of all eight states within the ground- 
level are extracted. For the strongest pump used in our 
experiment, over 95% of the F — 2 population piles up 
in the mp = +2 level. The reconstructed population 
distribution is shown in the inset of each panel in Fig. [2] 

The detailed population measurements provide a sub- 
stantial dataset, which can be used to benchmark a the- 
oretical multilevel model (28[. In Fig. [3l we show a 
comparison between the experimental measurements and 
the results of our multilevel numerical model, for a wide 
range of pump powers. The fitted Voigt parameters of 
the measured spectra (T and a) were used as an input 
to the model. A single model parameter that determines 
the ground-level population equilibration was fitted to 
the data, cffectivly describing the diffusion-induced decay 
and wall collisions (decay due to spin-exchange collisions 
is negligible in our conditions). It is evident that our 
model describes well the population distribution within 
the ground levels for five decades of pump power. 

When the pump power is very weak a thermal distribu- 
tion within the ground-state is observed. As the pump 
power is increased, the medium is polarized and more 
population accumulates in the higher mp levels of the 
F = 1 manifold (triangles and squares in Fig. [3]). At 
higher pump powers, the total population of the F = 1 
manifold is depleted due to off-resonant pumping. All 
these phenomena are well described by the multilevel 
model with a single fit parameter. Remaining discrepan- 
cies result from magnetic field inhomogeneity and from 
residual 85 Rb isotope in the cell (about 2%). It is in- 
teresting to note that an extremely efficient pumping to 
the maximally polarized state, \F = 2,mp — 2), can be 
achieved by tuning the frequency of the pumping laser 
in between the ground-level hyperfine transitions. We 
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verified this experimentally by locking the laser to the 
F = 3 — > F' = 2 transition of 85 Rb, and found that more 
than 90% of the ground-level population was pumped to 
the maximally polarized state for a pump intensity of 
about 100 /iW/mm 2 . 



B. Experiment II: Optical pumping with n 
polarized light 

We now turn to implement our reconstruction method 
for the case of CPT-based atomic frequency reference. 
One of the well known limitations of these frequency ref- 
erences is optical pumping of the atomic population from 
the clock-states \F = 1, mjr = 0) and \F = 2, mp = 0) 
Q. An approach previously investigated to resolve this 
problem is the repopulation of the clock states using 
additional 7r-polarized light field 0, [24| . To study this 
repopulation technique, we pump the ensemble with a 
RF-modulated 7r-polarized field, as depicted in Fig. [Tp. 
We use a current-modulated VCSEL at a modulation fre- 
quency of 3.0 GHz. The laser frequency is tuned such 
that the +1 and -1 sidebands are resonant with the 
F = 1 ^ f ' = 1 and F = 2 ^ f ' = 2 transitions, 
respectively. The laser is linearly polarized with a PBS 
and directed into a vapor cell containing isotopically pure 
87 Rb and 10 Torr of nitrogen buffer gas, at temperature 
of 58°C. The vapor cell is magnetically shielded and a 
constant magnetic field of about lOOmG is applied to set 
the quantization axis along the pump beam polarization 
direction (to induce it transitions). The probe beam, 
used for the population distribution measurement, is ap- 
plied along the magnetic field axis. 

Fig. |4] presents the reconstructed populations as well 
as the calculated results of our multilevel model. The 
\F = l,mtiF = 0) clock-state accumulates excess popula- 
tion as the pump power is increased, reaching a maximal 
population of twice the thermal equilibrium population. 
The populations in the mj? = ±1 states show symmet- 
ric behavior, demonstrating that a pure 7r-pumping was 
applied to the sample. As the pump power is increased, 
the population of the F = 2 manifold decreases, since 
the F = 2 — » F' transitions are more efficient than the 
F = 1 — > F' transitions. To emphasize the destruc- 
tive effect of the non-resonant excitations (blue dashed 
arrows in Fig. QJ), we have performed an additional cal- 
culation, in which the non-resonant excitations were ex- 
cluded (dashed lines in Fig. [4j. It is evident that the 
non-resonant excitations pose a limit on the efficiency of 
the 7r re-pumping technique. 



CONCLUSIONS 
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FIG. 4: (color online) The reconstructed populations ver- 
sus the pump intensity in experiment II, showing the agree- 
ment between the measured population (symbols), a numer- 
ical model (solid lines) and a numerical model that neglects 
non-resonant transitions (dashed lines) 



netic fields are weak, and the degeneracy created by the 
optical-transitions width is lifted by utilizing the proper- 
ties of the dipole matrix elements. As an intuitive quan- 
tity characterizing the state of the ensemble, the popula- 
tion distribution is most natural in monitoring dynamic 
processes, and in analyzing the initial quantum state in 
precise experiments [25j |. For the Dl transition of 87 Rb, 
the scheme reconstructs the detailed population distri- 
bution of the F = 1 manifold and the total population 
of the F = 2 manifold. In the framework of linear spec- 
troscopic measurements, the population distribution of 
the F = 2 level is fully reconstructed only if one can re- 
late the total population of the F = 2 manifold to one 
of the maximally polarized states. In order to recon- 
struct the full density matrix, a tomographic scheme (as 
in Q) could be applied. For that, the population distri- 
bution measurements, as presented in this work, should 
be carried out along AF + 1 different directions. We have 
used the scheme in two optical-pumping experiments and 
demonstrated its strength in describing various optical 
processes, such as resonant and non-resonant excitations 
and depolarization, and in verifying and calibrating the 
parameters of an elaborated numerical model. We fur- 
ther employed the calibrated model to calculate the im- 
provement of the contrast of a CPT based frequency ref- 
erence. 



We have presented an all-optical scheme to measure 
the populations within the ground-state of multilevel 
atoms and have demonstrated its implementation for a 
87 Rb vapor. In our scheme, the required light and mag- 
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The two eigenvalues remain small also if absorption mea- 
surements of a 7r-polarized probe are added to the recon- 
struction procedure. 

Our numerical model describes the dynamics and the 
steady-state solutions of the atomic density-matrix, 
within a Master-equation formalism, and was previously 
used in [15], H(| . The model incorporates all relevant sub- 
levels, e.g. the 16 sublevels of the rubidium-87 Dl transi- 
tion, and resonant coupling fields. Decay and decoherence 
processes, within the ground and the excited manifolds 
and between them, are effectively described by a suitable 
Lindblad form. We use 11 velocity groups to model the 
doppler broadening. 



